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Construction algorithms for plane nets in base b 
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Abstract 

The class of (0, m, s)-nets in base b has been introduced by Niederreiter as 
examples of point sets in the s-dimensional unit cube with excellent uniform dis¬ 
tribution properties. In particular such nets have been proved to have very low 
discrepancy. This property is essential for the use of nets in quasi-Monte Carlo 
rules for numerical integration. In this short note we propose two algorithms for 
the construction of plane (0, m, 2)-nets in base b. 
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1 Introduction 


In many applications, notably numerical integration based on quasi-Monte Carlo rules, 
one requires very uniformly distributed point sets in the unit-cube. The quality of the 
distribution is usually measured by the discrepancy. For a point set V = {a^i ,... ,X]\f} 
in [0,1)* and subintervals J C [0,1)^ the local discrepancy is dehned as 


A(J;P) 


£ {1,..., iV} : Xn & J} ^ Yohji 
N ^ ’ 


The discrepancy of V is then dehned as 


Dw(F)=sup|A(J;F)|, 

j 

where the supremum is extended over all intervals J C [0,1)^. Point sets with discrep¬ 
ancy of order Djq <C (logiV)^“^/Al are called “low-discrepancy point sets”. 

One class of point sets with excellent distribution properties are (f, m, s)-nets in base 
b as introduced by Niederreiter [5] (see also my- 

Let s,b & N, b > 2. An elementary 6-adic interval (or box) is an interval of the form 

TT hr G 
j=q ^ 

with hi,..., ds e Mo and Oj G {0,1,..., b^^ — 1}. 
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1 






Definition 1 (Niederreiter) Let b,s,m,t be integers such that b>2, s>l,m>0 
and 0 < t < m. A 6™-element point set V in the s-dimensional unit cube is called a 
{t, m, s)-net in base b if every elementary 6-adic interval of volume 5*“”* contains exactly 
b^ points of V. 


Thus a {t, m, s)-net in base 6 is a fe^-element point set V in the unit-cube for which 
the local discrepancy satisfies A{E]V) = 0 for all elementary 6-adic intervals E of 
volume The smaller t is (in the optimal case it is 0), the more demanding is this 

condition. This property for the local discrepancy is transferred in some sense also to 
arbitrary sub-intervals of [0,1)® which still have very low local discrepancy, although it 
cannot be zero in general. This is reflected in the discrepancy bounds which are of the 
form 


Dbm{V) <^s,b ^ 6 ^ ^ 


Um 


N 


if AT = 6" 


and if P is a {t, m, s)-net in base b (see one of [21IH El E])- 


Motivated by the discrepancy bounds the parameter t in the definition of nets is 
called the “quality parameter” and, as already indicated, a quality parameter which is 
as small as possible would be appreciated. Unfortunately, the optimal value t = 0 is 
not achievable for all possible choices of parameter pairs {b,s). It is well known that a 
(0, m, s)-net in base b can only exist if s < 6 -|- 1. This follows easily from the following 
proposition: 

Proposition 1 There is no (0, 2,b + 2)-net in base b. 

This result is very well known and there are several proofs available in the literature 
(see, e.g., [21 111 El E])- Nevertheless we give here a short and new proof which is based 
on arguments from Graph Theory and which might bring some new aspects into the 
theory of {t, m, s)-nets. 

Proof. Suppose to the contrary that there is a (0, 2,b + 2)-net in base b. Then its 
projection to the hyperplane orthogonal to the (6-|-2)nd coordinate axis forms a (0, 2, b+ 
l)-net. To every point in this net we may assign one of the colors {0,..., 6 — 1}, namely 
the first digit of the {b -|- 2)nd coordinate of the point in the original net. 

Now we construct a graph G = (U, E) by setting V equal to the points in the 
projection (hence |U| = 6^) and letting {^ 1 ,^ 2 } G U iff ui and V 2 lie in the same 6-1-1- 
dimensional elementary box of volume b~^. According to our construction two adjacent 
vertices have different colors. 

Every point in V is contained in 6-1-1 such elementary intervals and every such interval 
contains exactly 6 points. Thus every vertex of G has degree r = (6-1-1) (6 — 1) =6^ — 1 = 
|U| — 1. Thus G is isomorphic to the complete graph with 6^ vertices, for which we have 
found a coloring by 6 colors. But this is of course impossible. □ 

Usually, (f, m, s)-nets are constructed with the digital method which is based on m x 
m-matrices over a finite commutative ring with 6 elements (one matrix per coordinate). 
If 6 is a prime power, then there always exists a (0,m,6 -|- l)-net in base 6. More 
information about {t, m, s)-nets can be found in [21 lH El E]- 
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In this short note we will mainly be concerned with the two-dimensional case. Here 
one example of a (0, m, 2)-net in base b is the well known Hammersley net in base b 




tm ^ tm—1 ^ 


62 




Combining resnlts by Dick and Kritzer [2] and by De Clerck [T] we obtain; for every 
(0, m, 2)-net V in base b we have 




1 

b^ 


Cbm + 9 


where Cb ; = 


6+1 

6-1 


if 6 is even, 
if 6 is odd. 


( 2 ) 


So the discrepancy of a (0,m, 2)-net is of order of magnitude 0{{\ogN)/N), where 
N = 6"* = \V\. According to a celebrated result by Schmidt [7j this order is the best 
possible for the discrepancy of any A^-element plane point set. 


Our aim is to present two interesting construction algorithms which can in principle 
construct every (0, m, 2)-net in base 6. This way we construct point sets with optimal 
order of star discrepancy. 


2 The first algorithm 


If we are given a (0, m, s)-net in base 6, say V = {aii, X 2 ,..., Xb^}, then each point Xn 
belongs to a 6-adic 6"^ x 6”^ box of the form 


n 

i=i 


Uj{n) Uj{n) + 1' 


(3) 


where Uj{n) = [6™'x„jJ whenever Xnj is the jth component of Xn. It is elementary to see 
that the (0, m, s)-net property remains valid if one shifts the elements of a (0, m, s)-net 
within their corresponding 6-adic boxes of the form (151) . 

So, instead of constructing point sets, we are now going to construct sets of 6-adic 
boxes of the form (j5|). 

For a given box X = 11^=1 be the set of all 6-adic boxes of volume 

b~^ which contains X as a subset. 


Sm{X) = {E : E is a, 6-adic box of volume 6 with X C E}. 

We propose the following algorithm for the construction of hnite sequences of 6-adic 
boxes X = ^): 
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Algorithm 1 Construction of a finite sequence of boxes of the form ([2]) 
1: Input: base b, resolution m 
2: Set n = 1; 

3: Set Ui = [0,1)2; 

4; repeat 

5: Choose an arbitrary box 

A„ = n 

i=i 

6: Set n = n + 1; 

7: Set 

Un = Un-l \ U 

EGSmi^n—l) 

8; until Un = ^ 

9; return Xi,X 2 ,...,Xn 


Uj{n) Uj{n) + 1 \ 


An example of a construction according to Algorithm [T] is illustrated in Figure [TJ 

Theorem 1 Algorithm [I] terminates after exactly 6 ™ steps and the outcome yields a 
(0, m, 2)-net in base b. In particular, the so constructed point set satisfies the discrepancy 
bound ([ 2 ]). 

We split the proof into several short lemmas: 

Lemma 1 Algorithm]^ terminates after at most 6 ™ steps, i.e., n ^b"^. 

Proof. Each Xj belongs to exactly one fe-adic interval of the form 

[O’l) forae{0 ,l,...,6--l}. 

Hence the result follows since 


b^-i 


u 

a=0 


a a + 1 

5m’ Jjm 


X [0,1) 


[ 0 , 1 ) 2 . 


□ 


Lemma 2 If Algorithm [I] terminates after exactly steps, then the output seguence 
Xi, X 2 ,..., Xbm constitutes a (0, m, 2)-net in base b. 

For the proof we need the following easy lemma which we state for arbitrary dimen¬ 
sion s: 

Lemma 3 1. There are exactly b^ many s-dimensional b-adic intervals with 

volume 6 “™'. 

2. Every b-adic interval of the form X = 11^=1 [^ 5 -^^^) is contained in exactly 
s-dimensional b-adic intervals with volume b~'^, i.e., \8m{X)\ = 


4 















Figure 1: Algorithm [T] for b = s = 2 and m = 3; The chosen Xj,... ,Xs are colored 
black. The corresponding 6-adic boxes are in dark gray. 
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Proof. 1. The requested number is given by 


oo Pl-1 

E E---E 1 

di,...,<Zs=o ai=0 as=0 

f^iH- \-ds=m 


_ ^ -hrfs _ 


'm + s — V 


di ,...,ds=0 

di-\ - \-ds=rr] 


m 


2. For fixed (di,..., dg) G Nq with di + ■ • ■ + = m there is exactly one choice 


(oi, ... ,as) e JJ{0, 1 ,..., 6'^^ — 1} 

j=0 


such that 


n 

i=i 


Uj Uj + 1 


6™’ 6™ 


t=i 


cij cij T 1 


Jjdj ’ ^dj 


Hence the requested number is exactly the number of (di,..., d*) G Mg with di + 
■ ■ ■ + dg = m and this is 


□ 


Now we give the proof of Lemma El 

Proof. Assume that Algorithm!!] constructs the finite sequence Xi, X 2 ,..., Xbm. Accord¬ 
ing to Lemma El we have \8m{Xi)\ = = m -|- 1. Each Aadic interval E G EmiXf) 

contains exactly one element of Xi,... namely X^, which is the correct portion. 

Of course, 

8m{Xt) n 8m{Xk) = 0 whenever £ ^ k. 

Hence the number of Aadic boxes of volume 6“”* which contain exactly one of Xi, X 2 ,..., Xf,m 
is given by 

J2\8m{X,)\ = b^{m + l). 

£=l 

This is already the number of all 6-adic boxes in dimension 2 with volume b~^. Hence, 
every 6-adic box in dimension 2 with volume 6“"* contains exactly one of Xi, X 2 ,..., X^m. 
This means that Xi, X 2 ,..., Xf,m constitute a (0, m, 2)-net in base b. □ 

Lemma 4 AlgorithmUl terminates after exactly 6™ steps. 


Proof. Assume that Algorithm [T] terminates after the step n < 6™. Then there exists a 
A; G {0,1,..., 6™ — 1} such that 




k k + 1 ' 


6™’ b^ 


X [0,1) for all £ = 1, 2,..., n. 


We may assume without loss of generality that k = 0. Hence Xi ^ [0, x [0,1) for 
all £ = 1, 2,..., n. 

Now consider elementary boxes of the form 


0, 


<m—l 


X 


k k + l' 
V b 


for /c = 0,1,..., fe — 1. 


(4) 
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Assume that all of these b boxes contain one of Xi,..., Then 


u 


h=l 


h h + l\ 

5m ’ Ijm j 



( 5 ) 


contains b intervals Xj. According to the pigeonhole principle there must be an interval 
among the 6—1 intervals in union ([S]) which contains two Xj's. This however is impossible 
dne to the dehnition of the algorithm. Thus we have shown that at least one of the 
intervals in (jTj) does not contain a Xj. 

Again it is no loss of generality if we assnme that Xj ^ [0, x [0, ■^) for all 

j = 

Now we consider elementary boxes of the form 


0 , 


1 

'jjm—2 


X 


k k + l\ 

62 ’ 62 ) 


for fc = 0,1, ■ ■ ■ 5 6 — 1. 


( 6 ) 


Assume that all of these 6 boxes contain one of Xi,..., X^. Then 


b-l 


u 

h=l 


h 6, + 1 \ 

_ _ I ^ 

1 ’ 'jjm —1 j 



(7) 


contains 6 intervals Xj. According to the pigeonhole principle there must be an interval 
among the 6—1 intervals in union ([7]) which contains two X^-’s. This however is impossible 
due to the dehnition of the algorithm. Thus we have shown that at least one of the 
intervals in (jB]) does not contain a Xj. 

Again it is no loss of generality if we assume that Xj ^ [0, x [0, for all 
j = l,...,n. 

If we continue this process we hnally hnd that Un contains a 6-adic box with side 
length b~^. For example if in each step k = 0 then this is the box [0, Therefore 
there is space to choose a further interval X„+i. □ 


3 Algorithm [T] for arbitrary dimension 

In principle we can formulate Algorithm [T] also in arbitrary dimension s >2: 
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Algorithm 1’ Construction of a finite sequence of boxes of the form ([3]) in dimension s 
1: Input: base b, number of dimensions s, resolution m 
2: Set n = 1; 

3: Set Ui = [0,1)*; 

4: repeat 

5: Choose an arbitrary box 

x„=n 

i=i 

6: Set n = n + 1] 

7: Set 

Un = Un-l \ U 

EG£rn{^n—l) 

8: until Un = ^ 

9: return Xi,X 2 ,...,Xn 


Uj{n) Uj{n) + ^ 

’ Ifn 


Then Algorithm[r]still terminates after at most 6™ steps (same proof as for Lemma[T]) 
and if it terminates after exactly 6™ steps, then the output again constitutes a (0, m, s)- 
net in base b (same proof as for Lemma [2]). The problem is that we cannot guarantee that 
the algorithm runs until the 6™th step, i.e., we do not have a counterpart of LemmalUfor 
dimension s > 3. Even for s = 3 and 6 = 2 there are instances where the algorithm stops 
before 6"*. An example is illustrated in Figure El A necessary condition in order that 
Algorithm [Tl runs until n = 6™ is of course that s < 6 + 1 since otherwise a (0, m, s)-net 
in base b cannot exist. 
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Figure 2: Situation in dimension s = 3 and base 6 = 2 where Algorithm [P] stops already 
after 2 steps. Each square is a slice of the cube of height 1/4. 


4 The second algorithm 

The second proposed algorithm is a recursive construction in two dimensions which 
mimics the construction of Hammersley. It relies on the observation, that every (0, m, 2)- 
net in base 6 with m > 2 induces a (0,m — l,2)-net on each of the 6 rectangles with 
exactly one side length equal to 6“\ that partition the unit square. 

The algorithm we propose synthesizes a (0, m, 2)-net by (recursively) generating 6 
(not necessarily different) (0, m — 1, 2)-nets, scaling and shifting them and adjusting the 
points in a way that guarantees that the net-property is satished. 

Let Q{b"^) := {0, ..., be the set of 6-adic rationals with denominator 

not larger than 6™. Recall that any 6-adic 6“™ x 6“™ box is uniquely described by 





































the coordinates of its left lower corner. So we assume in the following without loss of 
generality that the elements of a (0,m, 2)-net in base b belong to Q(6™') x Q(6™’). 

Let Sb denote the set of permutations of the b elements {0, ... ,b — 1}. 

We also need the following two mappings: 

• for 6 G N, 6 > 2, let A;, : —)■ R^ be the linear scaling by b~^ in the direction of 

the hrst coordinate, i.e., Ab{x,y) = {x/b,y), and 

• for m G N, and permutations tto, ... ,7rb^-i £ 'Sb let ipb,m '■ Q{b'^) x Q(6™'“^) —)■ 

X Q(fe™), 


'il^b,m{x,y) := {x,y + nbm-iy{xi)b where Xi = [bx\. 


Algorithm 2 Recursive construction of a (0, m, 2)-net in base b 
1: Input: base b, hnal resolution m 
2: SetPo = {(0,0)}; 

3: for n = 1 to m do 

4: Choose permutations ..., G Sb] 

5: Set Vn ■= Ab{Vn-i + (jjO)); 

6: Set Vn := i>b,n{'Pn) 

7: end for 
8 : return Vm 


Theorem 2 The output Vm of Algorithmic is a {0,m, s)-net in base b, and therefore 
its discrepancy satisfies the bound ([2]). 

The proof of this result follows from the observation that Vq is a (0, 0, 2)-net in base 
b and by the following proposition which is slightly more general than necessary. 

Proposition 2 Let b,m E N, b > 2. Let b (not necessarily different) (0,m — l,2)-nefs 
Vj in base b for j = 0,... ,b — 1 be given as well as 6"*“^ permutations ttq, ..., 7ibm-i^i G 
Sb- Construct a set V of b^ points as follows: put 

b-l 

V := [j Ab{Vj + (j, 0)) and V := ipb, 7 n{V)- 

j=0 

Then V is a (0, m, 2)-net in base b. 


Proof. We need to show that every elementary interval of volume b contains exactly 
one element of V. Consider hrst an elementary interval I of the form 


I 


QjI 

Oi + 1 \ 


b^ ) 


0-2 

Jjm—k ’ 


02 + 1\ 
^m—k J 


with oi G {0,..., 6^ — 1}, 02 G {0,..., fe™ ^ — 1} and with /c G {1,..., m}. We may thus 
write 


/ — Ab{J + (h 0)), 
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where i = [aih and where J is the elementary interval 

02 a2 + l\ 

ym—k ’ Ijm—k J 

of volume This interval J contains precisely one element of Vi and hence / 

contains precisely one element of Ai,{Vi + (*,0)). On the other hand we have Ai^[Vj + 
(j, 0)) n / = 0 for j 7 ^ i. Thus I contains precisely one element of V. Since the second 
coordinate of this element belongs to and since - 06 ,m adds at most ^ to the 

second coordinate, I contains precisely one element of V. 


J = 


ai 


— ib^ 


-1 


ai 


— ^ + l' 


bk 


-1 


bk 


-1 


X 


It remains to consider the case that I is an elementary interval of the form 

0-2 1^2 + 1 


/ = [0,1) X 


Let j = [a 2 b~^\. By construction, each of the b elementary intervals 


+ 1 


6 ’ b 


X 


j J + 1 


Jfn—l ’ Jjm—l J ’ 


for £ = 0,..., 6 — 1, contains precisely one element of V. But by construction of 
precisely for one of these elements the second coordinate equals 026 “”*. □ 


Remark 1 If we choose the identity for every permutation in Algorithm El then we 
obtain the Hammersley net in base b as output, i.e., Vm = 'Hmfi (see (fT|)). 
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